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Abstract 

We investigate the effect of event-by-event fluctuations of the fireball created in high energy nuclear collisions 
on hard probe observables. We show that spatial inhomogeneities lead to changes in the nuclear suppression factor 
of high momentum hadrons which can be absorbed in the quenching strength q. This can increase the theoretical 
uncertainty on extracted values of q by up to 50%. We also investigate effects on azimuthal asymmetries V2 and 
dihadron correlation functions. The latter show a promising residual signal of event-by-event quenching that might 
allow us to estimate the size of spatial inhomogeneities in the fireball from experimental data. 
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High momentum quarks and gluons are a convenient probe of the quark gluon plasma bubble formed in high 
energy nuclear collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). Those 
hard partons are produced in rare hard processes in nuclear collisions. The process how such quarks and gluons 
convert into jets of hadrons in the vacuum, e.g. if the they are produced in p + p collisions, is fairly well documented. 
On the other hand we are still working to fully understanding how the same process plays out in a medium whith 
which those partons interact. Over the years several ways to compute the energy loss of energetic partons in quark 
gluon plasma have emerged E|2][3]|4), using different approximations to master the full complexity of the problem 
(see [5,6] for recent reviews). 

Calculations based on these models fare well with data from RHIC on nuclear modification factors Raa for leading 
hadrons if the overall quenching strength can be adjusted as a parameter. This parameter is often chosen to be q, 
the average momentum square transferred to the high momentum parton per mean free path. However, there are 
large discrepancies between the extracted values of q even if other details of the calculation, like the simulation of 
the space-time evolution of the quark gluon bubble are kept identical Q. Moreover, there is reason to believe that 
phenomenological details in the calculation, e.g. from fireball evolution or early-time behavior have a significant effect 
on extracted values of q. 

We will argue here that inhomogeneities and fluctuations in the transverse spatial distribution of both the medium 
(the quark gluon plasma) and the hard processes have an effect on q as well. On the other hand, we will show that 
precision measurements of hard probes can provide valuable experimental constraints on such inhomogeneities [8 1. 

Let us consider a fast parton emerging from a point r in the transverse plane with zero longitudinal momentum, 
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Figure 1: Density of nucleon-nucleon collisions in the transverse plane in Au+Au collisions with average impact parameter b = 3.2 fm generated 
with GLISSANDO. Left panel: A random single event. Right panel: Average over 500,000 events. 



traveling at a given azimuthal angle ip in the transverse plane. For a wide class of jet quenching models the average 
energy loss for such a parton traveling through a medium with local density p is determined by the integral 

WM = f dTlfipP+Tfy) (1) 

along its trajectory where t is the time parameter. The parameter ft determines the length dependence and would be 
equal to 1 in order to describe the consequences of the finite formation time known as the LPM effect ||5]|6). The 
impact of energy loss on the parton spectrum as a function of momentum p T and angle t/r is governed by the integral 
over all emission points r weighted with n(f), the probability density of jets emerging at point r, 

j d 2 r j dT/7t(>W+T^) (2) 

For the case of small energy loss Ap T <K p T the correction to the spectrum is directly proportional to this integral. 
For all other cases the relation might be more complicated, but it is still a monotonous function of the density product 
n(r)p(r + re^,). 

Now let us assume that the densities n(f) and p(r) fluctuate event-by-event around expectation values h(r) and 
p(r): n — h + Sn, p = p + dp. The effect on quenching from the density product, averaged over many events, can then 
be broken down into two terms 

{n(r)p(r + Te^)) = h(f)p(f + rfy) + R(r, f + re^) (3) 

where the first term corresponds to quenching in a smooth, averaged event, while the deviations from it are described 
by the correlation function 

R{A,A) = {Sn{A)Sp{m (4) 

In other words deviations of quenching when correctly calculated event-by-event and compared to averaged fireballs 
are sensitive to the correlation between the density of jet emissions and the density of the medium along the path of 
the jet. 

If we were able to measure these deviations we could gain access to the correlation function R, and hence estimate 
the average spatial inhomogeneities in the fireball, despite the fact that high transverse momentum (Pj) data at RHIC 
is averaged over many events. E.g. R would tell us about the average size of the typical granularity in events or the 
size of so-called hot spots. Fig.[T|shows the distribution n(r) represented by the number of nucleon-nucleon collisions 
for a typical Au+Au collision with small impact parameter {b « 3.2 fm) at RHIC energies together with the average 
distribution n(f). The distributions were calculated with the Glauber Monte Carlo GLISSANDO |9|. Fig.|2]shows the 
correlation function R(j) = R(r, r + re,/,) calculated with GLISSANDO events. One can clearly see a region of positive 
correlation around the emission point (small r) and a tail of anti-correlation extending to the boundary of the fireball. 

This behavior of R confirms the intuitive expectations. The integrals over t and r will pick up positive and negative 
contributions to the correlation function and it is hard to argue whether the net effect on the spectra is an enhancement 
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Figure 2: The correlation function tR(t) from a position f = 4 fm e y (f = corresponds to the center of the fireball) in the directions if/ = (solid 
line) and if/ = n/2 (dashed line) in central Au+Au collisions calculated with GLISSANDO. 

or suppression of quenching compared to the case of average fireballs. However, interestingly one can argue that 
the azimuthal asymmetry v'2 should decrease since the anti-correlation tails extend all the way to the boundary of the 
fireball 0. 

We have run numerical simulations using our PPM package |8| utilizing two different energy loss models: (i) a 
simple LPM-inspired model (sLPM) with AE = c sLPM /i(r, ijj) and (ii) the non-deterministic ASW/BDMPS model 0. 
In both cases we model the densities n and p using the binary nucleon collision densities from GLISSANDO. Note 
that the parameters c sLPM and c A sw can be interpreted as the local quenching strength per density, q/p. 

Fig.[3]shows results for Raa of neutral pions obtained with the ASW model for three different centralities compared 
to data from PHENIX iflOl . First casw is fitted using averaged events without inhomogeneities ("Averaged"), then 
the same value of casw is used in an event-by-event calculation ("Event-by-Event"). We clearly see a reduction in 
quenching for all centralities and momenta going from averaged events to realistic, fluctuating events. However, we 
can adjust c A sw from 1.6 GeV to 2.8 GeV to refit the Raa data on the same level of accuracy in the event-by-event case 
("Refitted"). Calculations using the sLPM energy loss model yield similar results [8]. We conclude that in absence 
of any knowledge about the spatial inhomogeneities in the fireball this introduces an additional uncertainty on the 
extracted value of q (or any version of this quantity normalized to a density) of at least 50%. 

The next interesting question is whether after refitting Raa any observable consequences of spatial inhomogeneities 
remain. We find that V2 at high P T is decreased going from averaged events to an event-by-event calculation with the 
same q, as predicted above. With the quenching strength used to refit Raa the results for V2 are still slightly below 
the results in averaged events, but the deviations (typically less than 20%) might be too subtle to be useful. On the 
other hand, triggered two particle correlations might be sensitive enough to put constraints on R. When going from 
averaged to event-by-event calculations using the same q the decreased suppression almost cancels between hadron 
pairs and trigger particles, leading to very similar curves. After readjusting Raa we find a rather large suppression of 
Iaa as shown in Fig. [3] 

It is also interesting to ponder the effects that inhomogeneities can have on azimuthal asymmetry coefficients other 
than \'2- We have analyzed a set of engineered events with given spatial anisotropies and have found that the Fourier 
coefficients v„ scale linearly with the (generalized) spatial eccentricities e„ in both the ASW and the sLPM energy 
loss models for n = 2, 3, 4. Using samples of GLISSANDO events we have indeed found Fourier coefficients up to 
n — 6. However except for V2 and V4 the magnitude of these coefficients is generally below 1%. Further details will 
be reported in a forthcoming publication |[T2l . 

Let us summarize. Inhomogeneities in the space-time structure of a quark gluon fireball have a potentially large 
effect on hard probes. One can define a correlation function R which encodes valuable information on the granularity 
and average magnitude of these fluctuations and which can potentially be extracted even in event-averaged measure- 
ments in heavy ion collisions. We have studied the function R within the GLISSANDO Glauber model and have 
investigated the effect of event-by-event quenching for two different energy loss models. We find that Raa increases 
and v>2 decreases event-by-event compared to averaged events if the same quenching strength q per density is used. 
However, we can not use this effect to constrain R since a simple redefinition of q can fit RHIC data quite well. 
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Figure 3: Left panel: Nuclear modification factor Raa for neutral pions in Au+Au collisions in three different centrality bins. Data from PHENIX 
1101 is compared to calculations using averaged events with q fitted (dotted lines), event-by-event calculations with the same q (dashed lines), and 
event-by-event calculations with q refitted to the data (solid lines). Right panel: Triggered correlation Iaa for central Au+Au collisions as 

a function of associated particle momentum Pj . The trigger particle momentum lies between 7 and 9 GeV. Data from the PHENIX Collaboration 
1111 is compared to our calculations for the averaged event (dotted line), the event-by-event calculations with unchanged q (solid line) and the 
event-by-event calculation with q refitted to describe Raa ■ 



However, in turn this means that there is a roughly 50% additional uncertainty on the extracted values of hatq which 
comes from our lack of knowledge about R. On the other hand we observe that two particle correlations like Iaa carry 
residual signatures of the correlation function R. If all other uncertainties were under control they could be used to 
experimentally constrain R and gain tomographic insight into the spatial structure of the fireball. We have also found 
a linear scaling of v n with e„ at large momentum and we find higher order Fourier coefficients in event-by-event jet 
quenching calculations which might be accessible in future measurements. 
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